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1  Introduction 


The  kinetic  Doi-Hess  theory  [6,  15]  has  been  a  useful  tool  for  modeling 
nematic  liquid  crystal  polymers  (TCP).  The  Doi-Hess  theory  uses  rigid  rods  to 
represent  the  nematogenic  molecules  and  describe  the  ensemble  with  an  orien¬ 
tational  probability  density  function  (pdf).  The  rotational  transport  equation 
for  the  orientational  pdf  is  given  by  a  nonlinear  Smoluchowski  (Fokker-Planck) 
equation.  The  Smoluchowski  equation  has  attracted  a  lot  of  attentions  from 
the  mathematical  society  [2,  3,  4,  7,  8,  9,  10,  11,  12,  13,  14,  18,  19,  23,  21,  17, 
24,  26,  22,  25]. 

In  [24]  we  showed  that  for  large  values  of  the  nematic  strength  and  small 
perturbation  there  is  a  solution  near  the  corresponding  unperturbed  pure  ne¬ 
matic  solution.  The  existence  was  established  using  complicated  free  energy 
arguments.  Furthermore,  only  existence  was  established  while  the  uniqueness 
was  not  shown.  Here  in  this  study  we  use  the  asymptotic  analysis  to  show  that 
for  large  values  of  the  nematic  strength,  the  Jacobian  is  nonzero.  Therefore, 
we  establish  the  extendability  (both  existence  and  uniqueness)  of  the  equilibria 
of  nematic  polymers  in  the  presence  of  small  perturbations. 


2  Nonlinear  systems  and  equilibrium  solutions 


We  briefly  present  the  mathematical  model  system  of  the  Doi-Hess  kinetic 
theory  for  homogeneous  flows  of  rigid  rodlike  nematogenic  molecules  immersed 
in  a  viscous  solvent  [6,  15].  By  “homogeneous”  it  is  meant  that  the  LCP 
orientational  distribution  is  uniform  in  space.  We  denote  the  orientational 
direction  of  each  polymer  rod  by  a  unit  vector  m.  For  pure  nematic  polymers, 
the  total  potential  consists  of  only  the  Maier-Saupe  interaction  potential 


t/(m)  =  — 6(mm)  ;  mm 


(1) 
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where  the  tensor  product  mm  and  tensor  double  contraction  A  ;  B  are  dehned 
as 


mm 


/  rriimi  mim2  mim^  \ 
m2mi  m2m2  rn2m^  , 
y  rrizmi  m^m2  mzmz  y 


A;B  =  Y^aijhiy  (2) 

In  (1)  b  =  3 A/2  and  A  is  the  normalized  polymer  concentration  describing 
the  strength  of  inter-molecular  interactions,  and  (mm)  is  the  second  moment 
of  the  orientation  distribution: 


(mm)  =  /  mm  p(m)  dm,  (3) 

d||m||  =  l 

where  p(m,  t)  is  the  orientational  probability  density  function  of  the  ensemble, 
i.e.,  the  probability  density  that  a  polymer  rod  has  direction  m  at  time  t. 
The  potential  (1)  has  been  normalized  with  respect  to  ksT  where  ks  is  the 
Boltzmann  constant  and  T  the  absolute  temperature. 

For  pure  nematic  polymers,  the  equilibrium  solutions  of  the  Smoluchowski 
equation  are  described  by  the  Boltzmann  distribution  [6]: 

1 

p(m)  =  -exp[-f/(m)],  (4) 


where  Z  =  exp[— f/(m)]  dm  is  the  partition  function  and  S  is  the  unit 

sphere. 

Let  us  select  the  coordinate  system  such  that  the  second  moment  (mm)  is 
diagonal: 


(mm)  = 


(  (m\)  0  ^  \ 

0  {nil)  0 

V  0  0  (mi)  y 

As  a  consequence,  the  Maier-Saupe  potential  can  be  written  as 


(5) 


f/(m)  =  —h{{ni\)ni\  -|-  {ml)ml  +  {ml)ml).  (6) 


The  most  signihcant  conclusion  for  pure  nematic  polymers  is  that  all  equilib¬ 
rium  solutions  are  axisymmetric  [7,  18,  23].  Since  not  all  equilibrium  solutions 
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of  a  perturbed  nematic  polymer  ensemble  are  axisymmetric,  to  facilitate  the 
discussion  of  the  extendability  ,  we  formulate  the  problem  without  using  the 
axisymmetry.  In  the  Boltzmann  form  (4),  an  equilibrium  solution  is  com¬ 
pletely  specihed  by  the  second  moment  (mm).  Because  of  the  constraint 
ml  +  ml  -|-  =  1,  an  equilibrium  solution  is  completely  specified  by 

Si  =  (mj),  S2  =  (ml).  (7) 

We  choose  the  coordinate  system  such  that  for  nematic  polymers  without 
perturbation  we  have  si  =  S2.  In  terms  of  (si,S2),  the  Maier-Saupe  potential 
can  be  expressed  as 

t/(m)  =  —b{sim‘l  +  S2ml  +  {l  —  Si  —  S2)ml) 

rl  3  ■ 

=  -b  2^S2  -  si){ml  -  mj)  +  {1  - -{si  +  S2))ml  +  const.  (8) 

From  (7),  the  nonlinear  system  for  (si,S2)  is 

Si  -  (mj)  =  0,  S2  -  (ml)  =  0,  (9) 

where  the  equilibrium  pdf  used  in  averaging  is 

p(m)  =  ^  exp  jfe  ^(s2  -  si){ml  -  mj)  +  (1  -  ^(si  -f  S2))ml  }  •  (10) 

The  form  of  the  equilibrium  pdf  (10)  suggests  us  to  introduce  a  new  pair  of 
unknown  variables 

[3  1  b 

Vi  =  b  l--(si  +  S2)  ,  V2  =  ^is2-si).  (11) 

Or,  equivalently,  the  nonlinear  system  for  {f]i,r]2)  is 

Fi{vuV2)  =  j-  - 1)  =  0, 

F2ivyV2)  =  j-^{ml-ml)  =  0,  (12) 

where  the  equilibrium  pdf  used  in  averaging  is 


p(m,  pi,p2) 


exp  [ri2{ml  —  ml)  +  pimg] 
exp  [772  (^^2  ~  ^1)  +  hi^i]  dm 


(13) 
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3  Jacobian  of  the  nonlinear  system 


Recall  that  we  have  selected  our  coordinate  system  such  that  for  pure 
nematic  polymers  we  have  772  =  0.  Let  us  calculate  the  Jacobian  of  system 
(12)  at  an  equilibrium  of  pure  nematic  polymer  where  772  =  0. 

We  start  by  hnding  the  partial  derivatives  of  pdf  (13): 


A 

dT]i 


p(m,  7/1,772) 


d  exp  [772(77^2  —  ml)  +  771 7n|] 
drji  Jg  exp  [772(7772  —  777 1)  +  7717773]  dm 


=  (777^  -  (r77^))p(m,7/i,r/2), 


(14) 


A 

dm 


p(m,  771,7/2) 


d  exp  [7/2(7772  —  777i)  +  7/1 7773] 

dm  Is  [^2(^2  ~  ''^1)  +  hi^i]  dm 


=  (7772  -  777i  -  (7772  -  777i))p(m,  7/1,  7/2).  (15) 

Then  the  partial  derivatives  of  (12)  at  an  equilibrium  (7/1  =  r,m  =  0)  are  given 
by 

■^Fi{m,V2)  =  {rnl)))  =  ^-^Yai{rnl), 

-^^Fi{r]i,m)  =  {rnl-rnl)))  =b,  (16) 

d  1 

^T2(pi,7/2)  =  {ml)))  =b, 

^T2(pi,7/2)  =  \-\{{ml-ml){ml-ml- {ml-ml)))  =  ^-]^{{ml-ml)‘^). 

In  the  above  calculations  all  averages  are  evaluated  using  the  pdf 

p(m,  Pi  =  r,  p2  =  0)  =  ^  exp(r7?73),  (17) 

where  r  satishes  the  equation 

(3  A)  -  1)  =  0.  (18) 

Now  we  write  (18)  in  a  more  explicit  form  using  spherical  coordinates.  We 
select  the  2:-axis  as  the  pole  of  the  spherical  coordinate  system.  The  pdf  (13) 
in  spherical  coordinates  is  of  the  form 

exp(r  cos^  0) 

277  Jq  exp(r  cos^  (/))  sin0d0 


p(</>, ») 


(19) 
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Substituting  this  expression  into  3{mf)  —  1  yields 

Jq{3  cos^  0  —  1)  exp(r  cos^  0)  sin  0  dcj)  fo(3u^  —  1)  exp(ru^)du 


3(7723}  -  1  = 


Jq  exp  (r  cos^  0)  sin  0  (i0 


Jq  exp(ru^)  du 


(20) 

where  a  change  of  variable  u  =  cos  0  is  applied.  Integrating  by  parts  yields 

2\  1  _  Id  exp(ru^)d(u^  "  ^)  _  o  "  u^)uexp(ru^)du 

Jq  exp[ru^)  du  Jq  exp{ru^)  du 

With  the  aid  of  (21),  the  equation  (18)  for  r  becomes 

1  Jq  u^(1  —  u^)  exp{ru^)du 


=  0. 


(22) 


_b  Jq  exp(ru^)  du 

Finally,  it  follows  from  (16)  that  the  Jacobian  of  the  nonlinear  system  (12) 
at  an  equilibrium  is  given  by 


d(FyF2) 


d(vi,V2) 


r/i=r 

V2=0 


13  f  2) 
-  -  -var(m3) 


(23) 


We  discuss  several  cases. 


Case  1:  Equation  (22)  has  a  trivial  solution  r  =  0,  corresponding  to  the 
isotropic  state.  In  the  isotropic  state. 


p(m)  =  {772‘^j)  =  (mj)  =  i,  =  ^,  i  ^  j, 


var 


(m3)  =  (m3)  -  (m3) 


(24) 


((m2  -  mi)  )  =  (m2)  -  2(mim2)  +  (m^  =  — 


5  15  15 

On  substituting  these  results  into  the  right-hand  side  of  (23),  we  obtain 


^3  12  15 

T  -  oVar(m3)  =  -  -  —  =  0  if  and  only  if  &  =  — , 


b  2 
1  1 


(25) 


^{{Tnl  -  Tnlf)  =  ^ 


—  =  0  if  and  only  if  5  =  — . 
15  2 


It  follows  immediately  that 

d{FyF2) 


diVl:V2) 


=  0  if  and  only  if  6  = 


15 


(26) 


ri\=r 

»?2=0 

15 


Therefore,  away  from  b  =  -y,  the  isotropic  state  can  be  extended. 
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•  Case  2:  For  r  7^  0,  equation  (22)  is  equivalent  to 


-  =  /(r),  where  /(r) 


Jq  u^(1  —  u^)  exp(ru^)  du 
/q  exp(ru^)  du 


(27) 


In  [23]  it  has  been  shown  that  the  function  /(r)  has  the  following  prop¬ 
erties: 


1-  /(0)  =  ^; 

2.  lim  f(r)  =  0  and  lim  f(r)  =  0; 

r— >+00  r— >— 00 

3.  /(r)  attains  its  maximum  at  r  =  r*  =  2.1782879748  where  f{r*)  = 

0.14855559992254  and  6*  =  ^  =  6.7314863965; 

f(r*)  ’ 

4.  For  r  <  r*,  f{r)  is  strictly  increasing;  For  r  >  r*,  f{r)  is  strictly 
decreasing. 

From  the  properties  of  the  function  /(r),  it  is  clear  that  equation  (27) 
has  a  solution  for  b  <  b*.  For  b  >  b*,  there  are  two  solutions:  rp{b)  >  r* 
and  Toib)  <  r*.  Here  we  use  the  subscript  “p”  to  refer  to  a  prolate  state 
and  the  subscript  “o”  to  refer  to  an  oblate  state  even  though  this  labeling 
is  not  completely  precise.  In  fact,  rp{b)  corresponds  to  a  prolate  state 
for  all  values  of  b  >  b*  while  ro{b)  corresponds  to  an  oblate  state  only 
for  b  >  7.5.  For  7.5  >  b  >  b*,  Toib)  corresponds  to  a  prolate  state.  As 
b  — >  -|-cx),  we  have 


lim  Tpib)  =  -|-cx),  lim  ro{b)  =  —00.  (28) 

6— >+00  6— >+00 

Now  we  look  at  the  asymptotic  behaviors  of  the  prolate  branch  (r  — 
-|-cx))  and  the  oblate  branch  (r  — >  — cx)).  For  mathematical  convenience, 
we  use  r  as  the  independent  variable. 

•  Case  2A:  The  prolate  branch  (r  — -|-cx)) 

Our  approach  is  to  expand  /(r)  for  r  — -|-cx).  First,  to  evaluate  the 
integrals  in  (27),  we  make  the  change  of  variables: 

V  =  1  —  u^,  u  =  Vl  —  V, 

dv  =  —2u  du,  du  =  — ,  dv. 

2^/T^ 
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Substituting  into  the  integrals  in  (27)  yields 

1 


/  exp(rn^)  du=  exp(r(l  —  v)) 
Jo  Jo 


2^1^ 


dv 


exp(r) 


rl  13 

/  exp(-rn)(l  + -n  + H - )  dv 

Jo  2  Q 


exp(r)  /111  3  1 

2  ^  1^0  ' 


/  n^(l  —  exp(rM^)  (in  =  /  (1  —  n)n  exp(r(l  —  n)) 
Jo  Jo 


exp(r)  /•!  /  I  3  2  \ 

J  exp(-rn)(l  -  n)n  [^1  + -n  + -n  H - J 

[  exp(— rn)  (v - - dv 

Jo  V  2  8  / 


2 

exp(r) 


2^/1^ 

■  ■  ■ )  dv 


dv 


exp(r)  f  1  1  3  1 

j^2  ^3  ^4 


+ 


By  using  these  asymptotic  results,  equation  (27)  can  now  be  written  as 


1  ,  ,  /o^  n^(l  —  n^)  exp(rM^)  dn  ^2 

b  Jq  exp{ru‘^)  du 


-  1  3  1 

—  (1 - + 


^  113  1 

1  +  —  ■  — 

2  r  4 


=  -(l 

r 


3  1  3  1 


+  ■■■ 


2  r  4 

The  order  parameter  is  given  by 

1  .  .  2\  \  ^  3  13  1 

2(3("i3)  = 

Now  we  calculate  var(m|)  and  ((m^  —  m\y).  We  note  that 
var(m3)  =  var(l  —  m\)  =  ((1  —  m\Y)  —  (1  —  rn'l)'^. 
Introducing  change  of  variable  u  =  cos  0,  we  get 

(1  -  ml)  = 


2-  Jq  sin^  0  exp(r  cos^  0)  sin  0  <70  /g^  (1  —  n^)  exp(rn^)  du 


Jq  exp(r  cos^  0)  sin  0  <70 


Jq  exp(ru^)  du 


{{l-miy)  = 


2\2\  fo  ^  exp(r  cos^  0)  sin  0  d(f)  /g^(l  —  n^)^  exp(rM^)  du 


fy  exp  (r  cos^  0)  sin  0  d(f) 


Jq  exp(ruy  du 


(29) 


2  2',2\  |/o’^sin^0  exp(rcos^0)  sin0(i0  1 /g^(l  —  n^)^  exp(rn^)  du 

/J"exp(rcos2  0)  sin0d0  2  /g^  exp(rn2)  dn 
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Note  that 


dv 


fl  rl 

/  (1  —  exp(rM^)  (iw  =  /  u  exp(r(l  —  w)) — -j= - 

7o  h  2vl  —  V 

=  ^  /‘  exp(-™)(„  +  + . . .)  +  1  +  . 

2  Jo  2  2 


[  {1  —  exp{ru‘^)  du  =  [  exp(r(l  —  w))- 
Jo  Jo 


dv 


2y/T^ 

exp(r)  /•!  ,  ,  ,  2  1  s  X  ,  exp(r),2  3 

=  I  oM-rv)(y  +  -v^  +  ...)*  =  -^(;5  +  ;i  + 

Substitution  of  (30)  into  (29)  leads  to 


{I -ml)  = 


Jq{1  —  u'^)  exp(rM^)  du 
Jq  exp(ru^)  du 


1.1 

;^(1  +  “  +  ■■■) 

t;  n  ■ 
;(i  +  2';  +  '") 


l.x  11 

—  “(1  +  ■  “  S - )) 

r  2  r 


((i-^^r)  = 


2  ,  3  1 

/o^(l  —  u^)^  exp(rM^)  du  ^  2  '  r  ' 


/q  exp(ru^)  du 


1.  11 

“(1  +  9  ■  “  + 

r  2  r 


2  .  1 
“  +  r  + 


{{ml-mir)  =  i- 


1.1 

=  19(1  +  r  — ) 


Hence 


1/1  3  1 

1  /^(l  —  u'^Y  exp(rn^)  du  ^3^1  3~  2  '  r  ^ 


Jq  exp(ruY  du 


1/1  1  1 

“(1  +  9  ■  “  d 
r  2  r 


-  -  -var(m3)  =  -  -  -  ((1  -  m3)  )  -  (1  -  m3 


b  2 


b  2 


,2\2 


1.  3  1  31  3 

“  ~  2  '  r  ~  4^  ^  ~  2 


2  ,  1  ,  1  ,  11 
;3(1  +  “  +  '")“;:7(1  +  2  ■  r  ^ 


1/1  3  1  31  ,  3  1  1 

-  7(1  ~  o  ■  -  -  TT?  - )  ~  n'  17(1  +  7  - ) 


2  r  4 


2  r2 


•  •),  (30) 


■■■)• 


■) 

■) 


■■■) 
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1  19  1 

=  -(1  —  3  - - T  ■  0  large  r, 

rjrt  lyt  ^  ly 


1  1  //  2  2N2V  1/.  3  1  3  1 

1  15  1 

=  -(1  —  2  - - T  ■  “^  +  ■■■)>  0  for  large  r. 

ry  y  ^  yZi 


Consequently,  it  follows  from  (23)  that 

^(^1,^2) 


11/.  1 
2  V*  + 


(31) 


diVl:V2) 


>  0  for  large  r. 


r]l=r 

V2=0 


Therefore,  for  large  r  (large  b),  the  prolate  branch  is  extendable. 

Case  2B:  The  oblate  branch  (r  — — cx)) 

We  use  the  same  approach  as  before  where  the  key  step  is  to  expand  /(r) 
for  A  =  — r  — >  +CX).  To  do  so,  we  introduce  the  change  of  variables 


dv. 


V  =  u  ,  u  =  \/V 


dv  =  2u  du,  du  = 


2v^ 


Then  the  integrals  in  (27)  turn  into 

fl  rl  2 

/  exp(rn^)  dn  =  /  exp(— An)  dv 
Jo  Jo  2,^/v 

_  1  r(i)  _  y?  1 


! _ y 


2  TA  2 

/  {1  —  u^)  exp{ru^)  du  =  /  (1  —  n)n  exp(r(l  —  n)) 

Jo  Jo 


2VT^ 


dv 


=  /  (1  —  n)n  exp(— An) 
Jo 


dv 


2^ 

—  2  J  (v^  —  v^J  dv 


TT  / 1 


1  /r(j)  _  r(j) A  ^ _ 

2  1a3/2  A5/2  ]  2  1 2  (-r)3/2  4  {-rf/^ 
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With  the  help  of  these  results,  (27)  can  be  rewritten  as 


=  f{r)  = 


1 

1 

'1 

3 

1 

Jq  u^(1  —  u^)  exp(rn^)  du 

(-r) 

_2 

~  4  ' 

(-r)_ 

i-r 


Jq  exp(ru^)  du 
3  1 


4  (-r) 

The  order  parameter  is  given  by 

1.  ,  9v  ^  r  13  1 

ysK)  - 1)  =  -  = 

In  order  to  calculate  var(m|)  and  ((m^  —  m\y).  we  observe  that 

/  cos^  0  exp(r  cos^  0)  sin  0 

^  Jq  exp  (r  cos^  0)  sin  0  <70 

Jq  u*  exp(rn^)  du  fo  exp(rn) dv 
Jq  exp(ru^)  du 

_  ir(i)  ■  (-r)-i  _ 


(32) 


\/n  1 

2  ' 


y/a  1 

2  ■ 


3 

4^2  ’ 


and 

(ml)  = 


2^  Jq  cos^  0  exp(r  cos^  0)  sin  0  dcj) 
fj  exp  (r  cos^  0)  sin  0  dcj) 

Jq  u^  exp(ru^)  du  fo  exp(rn)  dv 
Jq  exp(ru^)  du 


y/a  1 

2  ■  yf^ 


— r  2 


ypa  1 

2  ■  yf^ 

So  we  have 


var(m3)  =  (m3)  -  (m3)  = 


1 

2r 


4^2 


2r ' 


1 


((m 


2 

2 


1  /q^(1  —  u^Y'  exp(rn^)  du 

2  Jq  exp(ru^)  du 
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TT 


1  2 
2 


+ 


TT 


2(-r) 


+ 


(33) 


2 

Collecting  all  the  above  results,  we  find 


-  -  -varl^mgj  = 


1  3  1 

2  ~  4 ' 7^ 


3  1 

2  '  ^ 


>  0  as  r 


2r  2r2 

1  1  1 
^  //  2  2\2\  ^ 

6 -  “■))=(=;:) 


-CX), 


13  1 
2  ~  41^ 


1  (1 


2  V  2  2(-r) 


+  ■■■ 


+  ■  ■  ■  <  0  as  r  — — CX). 


4  4(— r) 

Thus,  from  (23)  we  arrive  at 

d(FuF^] 


(34) 


<  0  as  r 


-CX. 


d{vuV2) 

In  summary,  the  oblate  braneh  is  extendable  as  r  ^  —oo. 


r}i=r 

ri2=0 


4  Conclusions 

We  have  applied  the  asymptotic  analysis  to  show  that  the  equilibrium  so¬ 
lutions  of  the  pure  nematic  polymers  are  extendable  for  large  values  of  b  where 
b  is  the  nematic  strength. 
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